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About the Use of Real Dirac Matrices in 2-dimensional Coupled Linear Optics 
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The Courant-Snyder theory for two-dimensional coupled linear optics is presented, based on the systematic 
use of the real representation of the Dirac matrices. Since any real 4 x 4-matrix can be expressed as a linear 
combination of these matrices, the presented Ansatz allows for a comprehensive and complete treatment of 
two-dim. linear coupling. A survey of symplectic transformations in two dimensions is presented. A subset 
of these transformations is shown to be identical to rotations and Lorentz boosts in Minkowski space-time. 
The transformation properties of the classical state vector are formulated and found to be analog to those 
of a Dirac spinor. The equations of motion for a relativistic charged particle - the Lorentz force equations - 
are shown to be isomorph to envelope equations of two-dimensional linear coupled optics. A universal and 
straightforward method to decouple two-dimensional harmonic oscillators with constant coefficients by sym- 
plectic transformations is presented, which is based on this isomorphism. The method yields the eigenvalues 
(i.e. tunes) and eigenvectors and can be applied to a one-turn transfer matrix or directly to the coefficient 
matrix of the linear differential equation. 

PACS numbers: 47.10.Df, 41.75.-i, 41.85.-p, 05.45.Xt, 03.30.+p, 03.65.Pm 

Keywords: Hamiltonian mechanics, particle beam focusing, coupled oscillators, Lorentz transformation, Dirac 
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INTRODUCTION 

Even though there is continuous interest in this field 
(see, for instance^—), the treatment of coupled linear 
optics in two (or more) degrees of freedom has not yet 
reached the same level of generality, transparency and 
conceptual clarity as provided by the Courant-Synder 
theory for one degree of freedom. 

This article is about coupled linear optics as required to 
describe, for instance, the motion of charged particles in 
accelerators and ion beam optics. Even though ion beam 
optics is in principle tree-dimensional, often symmetries 
can be used to reduce the problem to the treatment of 
two-dimensional systems. In accelerators like cyclotrons 
or synchrotrons, the beam circulates in the horizontal 
plane and the electric and magnetic fields are symmet- 
ric with respect to this so-called median plane. In this 
case vertical motion is neither coupled to the horizontal 
nor to the longitudinal motion. But the dispersion of the 
bending magnets couples horizontal and longitudinal mo- 
tion. In other configurations, the transversal degrees of 
freedom are coupled with each other by solenoid magnets, 
but not with the longitudinal motion - or the longitudinal 
degree of freedom is "hidden" as the beam is not bunched 
but continuous. Therefore we will treat coupling in two 
dimensions like most authors^—. In an accompanying 
paper the problem of transverse-longitudinal coupling by 
space charge forces in isochronous cyclotrons is treated 
in linear approximation^ 2 -. This special type of coupling 
does not allow to apply the method of Teng and Ed- 
wards without modifications 12 . Other authors like Qin 
and Davidson restrict themselves to special forms of the 
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Hamilton function so that their treatment lacks general- 
ity. 

The real Dirac matrices (RDMs) have been known for 
a long time as the Majorana representation, going back 
to a paper by Ettore Majoran a 13 ' 14 . In this article, we 
will use the RDMs in a very practical way that has - to 
the knowledge of the author - not yet attracted much at- 
tention. The RDMs form a basis of 4 x 4-matrices with 
remarkable properties also - and maybe especially - in the 
context of classical mechanics. The RDMs are essential 
ingredients for a formulation of a theory of symplectic 
flow in two dimensions. They enable to survey all possi- 
ble symplectic transformations in an elegant and straight 
manner. 

The use of real instead of the complex Dirac matri- 
ces has several reasons: First, linear coupled optics is a 
classical theory and the relevant terms in the Hamilton 
function are real. Second, the RDMs are a complete basis: 
Any real 4x4 matrix can be written as a linear combi- 
nation of the RDMs. This in fact guarantees generality 
and completeness of the presented theory of linear cou- 
pled motion in two dimensions. And finally - the RDMs 
are discriminable with respect to all relevant structural 
properties: Each matrix is either symplectic or antisym- 
plectic, either even or odd, two RDMs either commute 
or anticommute. The use of the RDMs supports another 
clear distinction - real or imaginary. 

Furthermore the introduction of RDMs into classical 
mechanics may provide new insights into the relation- 
ship of Hamiltonian mechanics, special relativity and 
the Dirac equation. It is known for a long time that 
the mathematical formalism of Twiss-parameters and 
Courant-Snyder invariants can be applied to quantum 
systems 1 -. But no attempt has yet been made to apply 
the tools of quantum mechanics in classical mechanics. 

Nevertheless we emphasize that a lot of the presented 
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formalism, i.e. practically all equations that do not re- 
fer to other matrices than 70, may be applied in arbi- 
trary dimensions, if 70 is extended correspondingly. This 
holds especially for the concept of the symplex, that we 
introduce to identify the components of the "force ma- 
trix" . We will demonstrate the significance of the RDMs 
for the treatment of coupled linear motion, in the con- 
text of transfer matrices and eigensy stems. This leads 
in a natural way to the two-dimensional generalization 
of the Courant-Snyder theory. In the second part we 
demonstrate how the RDMs can be used in the context 
of the second moments using Poisson brackets. We show 
that both, the Maxwell-equations and the Lorentz force 
equations, can be formulated casually in terms of the 
RDMs. This brings up an analogy which we call the 
"electromechanical equivalence" . Effectively we use the 
isomorphism of symplectic transformations with Lorentz 
boosts and rotations in Minkowski space to introduce a 
physical nomenclature of symplectic transformations in 
two dimensions. Based on this equivalence and on the 
distinction of even and odd matrices we finally develope a 
general algorithm that allows to determine the symplectic 
decoupling transformation. Furthermore the algorithm 
enables to compute the eigenvalues and the eigenvectors 
of stable two-dimensional symplectic systems. 



I. REAL DIRAC MATRICES IN COUPLED LINEAR 
OPTICS 

The position and direction of a charged particle within 
a beam is usually described by its coordinates and an- 
gles relative to the reference trajectory. In case of two 
transversal coordinates, the position of the particle is de- 
scribed by x, x' , y, y' ', where x is the horizontal and y the 
vertical ("axial") coordinate. The dashes represent the 
derivative with respect to the path length s along the ref- 
erence orbit. In case of transverse-longitudinal coupling, 
a common choice of the coordinates is x, x 1 ,1,8 where I 
is the longitudinal coordinate with respect to the bunch 
center and 5 = p ~ Po is the relative momentum deviation, 

Po _ ' 

where p is the momentum of a specific ion and po is the 
average momentum. 

Since the formalism is related to classical Hamiltonian 
mechanics, we prefer to write qi and pi for the dynamical 
variables and we collect them in a vector 



1p = {q 1 ,p 1 ,q2,p 2 y j 



We could as well change the ordering of the variables and 
use for example 



'0 = (0'i,?2,Pi,P2) 5 



(2) 



This may (of course) have no influence on the physical 
situation and its description, but it leads to a different 
ordering of the RDMs. Four elements can be ordered in 
4! = 24 different ways, but since we do not distinguish 
the indices, there are 4!/2 = 12 different permutations of 



the elements in ip and corresponding systems of RDMs 
as listed in Tab. [TTT1 

We identify the Dirac matrix 70 with the symplectic 
unit matrix, which is usually denoted by /, J or S. In 
case of Eq. [TJ the "time direction" 70 is 



To 



10 
-10 
1 
0-10 



(3) 



in case of Eq. [5J the form is 

/ ° 


m = _ 1 



V -100 




(4) 



70 is a skew-symmetric matrix that squares to the nega- 
tive unit matrix 7q = — 1 , while in the usual definition of 
the Dirac matrices one has 7q = 1. The other three ba- 
sic matrices 71 ... 73 are defined by the requirement that 
they must anticommute with 70 and with each other and 
that they square to the opposite sign of 79. The signa- 
ture of the metric tensor g^ v , hence, is (—1, +1, +1, +1) 
- instead of (+1, -1, -1, 



= Diag(-1, 1,1,1) 

_ 7n 7»*+7i/ Jit 
2 



(5) 



The other matrices that we use to form the symplectic 
basis are 



7i = 



72 



73 



0-100 
-10 
1 
10 

/ 1 
10 
10 

V 1 

-10 
10 
0-10 
1 



(6) 



(7) 



(8) 



(X) The remaining matrices are defined by 
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= 70717273; 


715 


= 1 




74 


= 7o 7i; 


77 


= 714 7o 7i = 


72 73 


75 


= 7072; 


78 


= 714 7o 72 = 


73 7i 


76 


= 7073; 


79 


= 714 7o 73 = 


7i72 


710 


= 7i4 7o = 


71 72 73 






711 


= 714 7i = 


7o 72 73 






712 


= 714 72 = 


7o 73 7i 






713 


= 714 73 = 


7o 7i 72 







(9) 
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TABLE I. Properties of the real 7-matrices. The type encod- 
ing is V for vectors, B for bi-vector, A for axial vector, P 
for pseudoscalar and S for scalar and refers to the symplectic 
transformation properties. The subscripts t (s) indicate time- 
(space-) like components of a four-vector, respectively. The 
column labeled y a gives the sign s that fulfills the following 
equation: 7 a jx ja = sj x . The column labeled "s/a" indi- 
cates, whether 7^ is a symplex or an antisymplex, respectively. 
Accordingsly, the column labeled "S/A" tells, if 7^ is symplec- 
tic or antisymplectic, respectively. The column labeled "e/o" 
indicates, whether the corresponding 7-matrix is even or odd. 
Even matrices are non-zero only in the block-diagonal com- 
ponents, odd matrices are zero in the block-diagonal compo- 
nents. Finally, the column labeled 70 denotes, if a basis exists 
such that 7a, appears as the "time direction", i.e. if a basis 
exists in which 7^ plays the role of 70. This is only the case 
for antisymmetric matrices, which are either both a symplex 
and symplectic or none of it. 



have the familiar form 

Pi 

or, in vector notation, 



8H 
dpi 
dH 



ip = 70 Vw, H , 



(13) 



(14) 



where the dot represents the time derivative. It is well- 
known that the Jacobian Matrix M of a canonical trans- 
formation is symplectic, i.e. thati^ 



M 70 M T = 70 . 
A matrix M is antisymplectic, if 

M 70 M T = -70 . 
Eq. [15] includes also 

M T 70 M = 70 

M T = -7oM- 1 7o 
M" 1 = -7oM T 7o. 



(15) 



(16) 



(17) 



The product of symplectic matrices is symplectic. But 
note: Also the product of two anti-symplectic matrices 
is symplectic (see also Tab. [I]). 



Note that the definition deviates from the conventions in 
particle physics, where the product 70 71 72 73 is usually 
labeled 75. The matrices are explicitly given in App. |5] 
If A is an arbitrary real- valued 4x4 matrix, then A can 
be written as a linear combination of RDMs 
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^2 a k Ik ■ 

k=0 



(10) 



where the RDM-coefficients au are given by the scalar 
product 

Ofe = A • 7^ = — - — Tr( - ), (11) 

and Tr(X) is the trace of X. The RDMs form a group and 
are the basis of a vector space. The associated algebra is 
the real Clifford algebra C7(3, 1)^2. 



A. The Hamiltonian 

Linear coupled optics is characterized by a Hamilto- 
nian function of the classical harmonic oscillator form 



H = \V AV, 



(12) 



where the superscript "T" denotes a transposed vector 
or matrix. A is a (generally time dependent) symmet- 
ric matrix. The Hamilton equations of motion (EQOM) 



B. The Force Matrix and the Definition of a Symplex 

From Eq.[T2land Eq.[14]one derives the following linear 
EQOM: 



ip = 70 A -0 = F ip . 



(18) 



where F is the (generally time-dependent) force matrix. 
The antisymmetry of 70 and the symmetry of A yield: 



= A T 7(f 

= -A 70 

= 70 70 A 70 

= 70F70 



(19) 



In the following we call a matrix F that fulfills EQ. [T9la 
symplex (not "simplex" ) . Symplices are sometimes called 
" infinitesimally symplectic" or "Hamiltonian" 19 , but the 
author prefers a unique and short name. A matrix F a 
that holds 



-7o Fa 70 



(20) 



is called an antisymplex. 70 itself is a symplex as it is an- 
tisymmetric and squares to the negative unit matrix. By 
definition the basic matrices 71 ... 73 are also symplices. 
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The most important property of symplices is the super- 
position principle. According to this principle the sum of 
two symplices is a symplex: 



(F!+F 2 ) T = Ff + FJ =7oFi 7 o+7oF 27 o 
= 70 (Fi + F 2 )7o • 



(21) 



The superposition principle includes scalability: A sym- 
plex multiplied by a scalar is still a symplex. Given that 
the product of two symplices Fi and F 2 is also a symplex, 
then one finds: 



(22) 



7o(F 1 F 2 ) 7o = (F!F 2 ) T 

= F^Ff 

= 7oF 2 7o7 Fi7o 

= -7oF 2 Fi7o 

Fi F'2 = — F2 Fx 

FiF 2 + F 2 Fi = 0. 



The product of two symplices is a symplex, if (and only 
if) the symplices anticommute. Since the four basic ma- 
trices 70 ■ ■ • 73 anticommute with each other, the six pos- 
sible bi-vectors 7^ 7 M are symplices. 

Since a symmetric n x n-matrix A is described by v 
parameters with 



v = n(n+ l)/2. 



(23) 



the force matrix F = 70 A must contain the same number 
of independent components. In case of n = 4 we expect 
v = 10 force components. These are the four basic ma- 
trices, and the mentioned six bi-vectors 7„ 7^ = 74 ... 79. 
Hence, a general force matrix in two dimensions has the 
form 



F = E/< 



ft Ik ■ 



(24) 



k=Q 



C. Symmetric Products and Projectors 

The symmetric product of a matrix Fi, which is either 
a symplex or an antisymplex, and a symplex F 2 is again 
a symplex: 

(F 1 F 2 F 1 ) T = Fj F2 Fj 

= (±7 Fi7o)(7oF 2 7o)(±7oFi7o) 
= 7o(FiF 2 Fi)7 . 

(25) 

Since all 7-matrices are either symplices or antisymplices, 
any expression of the form 
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^ a k Ik F 7 fc 



k=0 



with arbitrary coefficients a& is a symplex, if F is a sym- 
plex. Tab. Ushows the result of j a ^fxja for a = [0, 10, 14]. 



The result of the operation 

F±7 Q F 7a 1 ^ 
F a = 2 = 2 ^ ^ la lk la ' ' \ ' 

k=0 

is a projection. For a = 14 for instance one has 

9 3 

I E /fe (7ft +7i4 7fc7i4) = E fklk 
*=o ft=o (27) 

5 E /fe (7ft -7i4 7fc7i4) = E fklk, 

k=0 fe=4 

that is, 71 4 separates the "vector" -components from the 
"bi- vector" components. 



D. The Transfer Matrix 

The solution of Eq. [15] can be written by a symplectic 
transfer matrix M(i, to)- 

^(t)=M(t,to)V(*o). (28) 
If the force matrix is constant in time, then 

M(t,io)=exp(F(t-t )). (29) 
The time derivative of EQ. I2"51is: 



#) = M(Mo)V>(*o) 
= FM(t,t )V(«o) 



(30) 



so that the following differential equation holds for M: 

M = FM. (31) 

Comparison with EQ. IT51 shows that the column vectors 
of the transfer matrix M are solutions of EQ. [18] If n 
is the dimension of the complete transfer matrix can 
be obtained by integrating Eq. [T5] n times, using the n 
euclidean unit vectors as starting conditions il>(to). 

If the transfer matrix is expressed by a time-dependent 
matrix <I? according to 



M = exp$= £ 



(32) 



then if (and only if) <fr and 3? commute, i.e. if 

= (33) 

Eq. [32] can be written as: 



M = + * + + ^ 

= $ exp«fr 
= $M 



(34) 
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so that in this case one finds 



The Definition of Coupling 



F = * 

*(Mo) - J F(t)dt = Fr, 



(35) 



to 



with t = t — to so that the transfer matrix can be written 
as 



M(t) = exp (Ft) . 

If €> and <& would anticommute, i.e. if 
+ = 0, 



(36) 



(37) 



then the square - and any even power of - $ would be 
constant and therefore one would find: 



M = E TOT= sinh (*)* _1 * 



fc=0 



(38) 



We define the matrices M s and M c according to 



■jyj _ exp (#)-cxp (-#) _ M-M -1 _ 

tv /r - exp (^)+cxp (— #) _ M+M 
1V1 C _ 2 — 2 



75 — = sinh (#) 



cosh (<&) 



(39) 



so that Eq. [37] would yield 



M c = 

m s = M,r'$. 



(40) 



In Sec. II HI it will become clear that Eq. [37] has to be 

rejected in the case of stable focused systems. 

Beam transfer lines and circular accelerators are typ- 
ically composed of guiding elements that provide a con- 
stant force matrix for a certain time (or better: length). 
Examples are dipole- or quadrupole magnets, drifts and 
bends. In this case, it is possible to express the transfer 
matrix as a product of transfer-matrices for the individ- 
ual elements: 

M(t n ,t ) = M(t n ,t n _i)...M(i k ,t fc _ x )...M(ti,to) 
= exp (F n r n ) . . . exp (F x n) . 

( 41 ) 

The transfer matrix is symplectic, if $ is a symplex: 



M 7o M T 



l~i fe! 



fe=0 



7o ( E 

Js=0 



fe=0 



\k-l. 



Ei ho E i~r A 7o irTo 



A'! 



k=0 
fc=0 

exp (<&) exp (— <I>) 70 

7o 



Li k\ 

k=0 



7o 



The exponential of a symplex is symplectic. 



(42) 



Before starting an investigation on decoupling there 
should be a clear definition of coupling. A possible (and 
typical) definition refers to the structure of the force ma- 
trix. As we defined two degrees of freedom qi and q2, 
the obvious form of decoupling is a block-diagonal force 
matrix: 



(43) 



where A and B are 2 x 2-matrices. 

There are less obvious forms of decoupling - consider a 
constant force matrix. The second derivative of the state 
vector tp is: 




tp = F ip = F 2 ip . 



(44) 



If the square of a constant force matrix is (block-) diag- 
onal, then the system can be regarded as decoupled in 
second order. The solutions can be found separately for 
each degree of freedom. In this case the coupling fixes 
the relative phases between the different degrees of free- 
dom - but it does not determine the functional form of 
the solution. 

A force matrix of the Dirac-form is an example: 

F = E-f +P171 +P272 +P3 73 (A r\ 
F 2 = -(Ei-pl-pl-pl)l. 1 ] 

Even though the odd component pi 72 is not block- 
diagonal, the second order differential equation is decou- 
pled. Another interesting example is a constant force 
matrix X of the form 



-y^ 7o~ 72 — 76 — 77 

X 2 = -711 
X 4 = 1 



(46) 



Here one finds that neither X nor X 2 or X 3 are block- 
diagonal. Nevertheless the fourth time derivative of tp is 
"decoupled" . 

We refer to systems as decoupled in n-th order, if the 
n-th order EQOM have the form 



ip = B if) , 



(47) 



where B is block-diagonal. In case of time-dependent 
force matrices, the second derivative is 



(F4 



HF</> 



(48) 



Hence we consider systems with time-dependent force 
matrices to be decoupled to second order, if the expres- 
sion F + F 2 is block-diagonal. 
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F. Symplectic Symplices 

Tab. U lists all RDMs with their main properties. Each 
RDM ist either a symplex or an anti-symplex, symplectic 
or anti-symplectic. Four RDMs are both - symplectic and 
symplex. If F is a symplectic symplex (SYSY) , then the 
combination of EQ. 02] and EQ. Q1J] yields: 



F 7o F 7 



F 70 70 F 70 
F 2 = 1. 



-F 2 70 = 70 



(49) 



A symplex is symplectic, if (and only if) its square equals 
the negative unit matrix. But note that symplectic ma- 
trices are in principle not scalable, i.e. without unit. 
Symplices are scalable and may therefore have arbitrary 
units. If the symplex F appears in the form of Eq. [TBI 
then it has the natural unit of a frequency or wavenum- 
ber if the dot is interpreted as time or pathlength deriva- 
tive, respectively. In practical problems of accelerator 
physics, it will always be possible to find a typical length 
or time interval that can be used to redefine the dot- 
derivative in such a way that the force matrix is unit- 
less. In consequence this means that force matrices which 
fulfill F 2 = — SI 2 1 with a constant are equivalent to 
SYSYs: 



dr = fldt 



(50) 



In this respect the Dirac-operator (Eq. |45|) is a SYSY. 

From Eq. 09] and Eq. [29] it is quickly derived that the 
transfer matrix of a constant SYSY F is given by 



M 



1 cos t + F sin t . 



(51) 




Eq. [51] is known in the Courant-Synder theory of 1-dim. 
ion beam optics. The unitless parameter t is called phase 
advance. If t represents the phase advance for one turn, 
then it is called tune. In the 1-dim. theory, the matrix 
F has the form 



(52) 



where a, j3 and 7 are the so-called Twiss-Parameters 2 ^,. 
According to the concept of the Sec. lIE) a constant SYSY 
is decoupled in second order and Eq. [51] shows that this 
definition is meaningful. Besides this a SYSY also has 
equal eigenfrequencies - as will be shown in Sec. II Gl - 
and in consequence there is a common phase advance for 
both degrees of freedom. 

Let U be a time dependent symplectic transformation, 
then we obtain for the transformed "spinor" ip: 



(53) 



The time derivative of EQ. [TBI is 

V> = F^ + FV' 



(54) 



The formal difference between Eq. [S3] and Eq. [5j] is that 
U is symplectic and F is a symplex. If U is also a symplex 
(or if F is symplectic), then these equations through light 
on the structural equivalence between a (time) derivative 
and the multiplication with a SYSY. If the dynamics of a 
system is described by a SYSY, then the (time) derivative 
is itsself a symplectic transformation. One could say that 
SYSYs are operators, which are equivalent to derivatives. 



G. Eigenvalues and Eigenvectors 

Eigenvalues and eigenvectors play an important role, 
if the system has some kind of self-feedback. Circular 
accelerators like storage rings are a simple example for a 
system with self-feedback. Another example are systems 
with a constant force matrix, so-called "constant focus- 
ing channels" . If A is the diagonal matrix containing 
the eigenvalues of F and E is the matrix of columnwise 
eigenvectors, then 



FE = EA. 

If E can be reversed, then 

F = EAE" 1 



(55) 



(56) 



In the simplest case of a constant force matrix one finds 
that 

M(r,0) = exp(EAE- 1 r) 



(E A E 



-1 



fe! 
k\ 



E- 1 



where 



= E 

fe=0 
00 

= E e 

k=0 

= E cx P (At)E- 1 
= EAfrjE- 1 , 



A(t) = cxp (Ar) 



(57) 



(58) 



is the diagonal matrix of the eigenvalues of M. 

The trace of the product of an antisymmetric and 
a symmetric matrix is zero. Hence F has zero trace. 
As similarity-transformations preserve the trace and the 
eigenvalues, we find Tr(X) = Tr(F) = and hence the 
sum of all eigenvalues is zero. The eigenvalues are either 
real or (two) pairs of complex conjugate values. A sym- 
plectic transformation is said to be strongly stable, if all 
eigenvalues of M are distinct and lie on the unit circle in 
the complex planed. This means that for stable (oscil- 
latory) solutions the eigenvalues of F are two conjugate 
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pairs of imaginary values: 



A 



Eg. [551 yields: 



Diag(i wi, - 



E A 2 E- 1 

2 i 2 
,U, +U1. 2 



lb}t,l<J2,—lW2) 



(59) 



-El 



I ^7i 2 )E- 1 
E712E- 1 . 



(60) 



This shows that F is (isomorphic to) a SYSY, if w 2 = 
w 2 . SYSYs are degenerate. The absolute values of all 
eigenfrequencies of a SYSY are equal. 



H. The Form of the Transfer Matrix 

Eq.[36l allows specific functional forms /(t) for the ma- 
trix elements - depending on the dimensionality and the 
properties of the force matrix. The force matrix of a 
strongly stable system has eigenvalues that are grouped 
in two pairs of imaginary values - the eigenfrequencies. 
Ring-accelerators always have this property - but not in 
each section. Ion beam transport systems are usually 
composed of sections with separate elements like dipole-, 
quadrupole or sextupole magnets, drifts, solenoids, and 
so on. These elements are characterized by their indi- 
vidual transfer matrices and are not necessarily "stable" . 
Only the product of the transfer matrices of all elements 
in a ring-accelerator has to be stable. Eq. [36] is actually 
computed as a series: 



M = exp (F t) 



(Ft)* 



(61) 



In case of pure RDMs, i.e. F = 7fe with k € [0...9], 
the functional form of the elements of the corresponding 
transfer matrices are exponentials of r - if we include (hy- 
perbolic) sine- and cosine forms. In the case of arbitrary 
symplices, other forms are possible: A square matrix F is 
called nilpotent, if F q = for some positive integer q > 1. 
A simple example is the "force" matrix of a drift, which 
is given by (see Eg. 11331 in Sec. Ill Gl below): 



F 

M 



70+76 



fc! 



£ 

fc=0 

1 + Fr 



(62) 



Since the maximal (non-zero) power of a n x n matrix F is 
n— 1, polynomials up to third order as well as products of 
polynomials and exponentials are also possible solutions 
of Eq. iU 

In case of a constant force or in case of a transfer ma- 
trix for a complete turn in a stable circular accelerator, 



the diagonalized transfer matrix can be computed from 
Eq. [55] and Eq. [SHI We introduce the abbreviations 



S, 
S s 

A c 
where 



cos (wi t)+cos (wo r) / - \ /a \ 

^ — ' ; = cos (w t) cos (Aw t) 

sin (w r) cos (Aw r) 

sin (cji r) — sin (u2 r) 

2 — 

cos (cji r)— cos (w2 t) 

2 — 



sin (uji r)+sin (cJ2 t) 
2 



cos (w r) sin (Aw r) 
— sin (w t) sin (Aw r) , 



Ul+^2 



Aw = 



(63) 



(64) 



so that the diagonal matrx A (Eq. [57] and Eq. [58]) can be 
written as 

A = S c 1 - i S s 73 - i A s 74 - A c 712 , (65) 

and the one-turn-transfer matrix of strongly stable sys- 
tems is according to Eq. [57] given by: 

M = S c l - i£ s E7 3 E- 1 

— i A s E74 E _1 — A c E 7 12 E _1 (66) 



Eq. [55] is the generalization of the Twiss-matrix for 2- 
dimensional systemsSi. It shows that transfer matrices 
usually have - in contrast to symplices - a scalar compo- 
nent as well as components which are antisymplices. The 
second and the third term of Eq. [55J have the same form 
as a force matrix (though different eigenvalues) . The last 
term vanishes in case of a degenerate system with equal 
eigenfrequencies. In this case the one-turn transfer ma- 
trix has the form of Eq. [51] 

In order to split the transfer matrix into the compo- 
nents according to Eq. [55J we make use of the method 
of projections, modified for this purpose, as follows: The 
inverse of a symplectic (transfer-) matrix M has the same 
eigenvectors as M: 



M- 1 = EA-'frjE- 1 

= E A(— t) E _1 = E exp (—At) E _1 

so that using Eq. [T71 yields: 

i(M±M-!) = E A(T)± 2 A '- r> E- 1 
= i(MT7oM T 7o ). 



(67) 



(68) 



If we recall Eq. [311 then we can determine the RDM- 
coefficients of the transfer matrix and the matrices M s 
and M c according to Eq. [TU] and Eq. [TT] Let nrik be the 
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resulting coefficients, then insertion into Eg. l68l yields: 
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E 



i(M + 7oM T 7 o) = £ m klk 

k=0 

A(r)-A(-r) F 1 
2 

= -iE s E7 3 E- 1 -iA s E7 4 E- 1 

15 

M c = i(M- 7o M T 7 )= ^ m fc7fc 



(69) 



fe=10 



E A(r)+A(-r) g _i 

S c l- A c E 7l2 E- 1 



A decoupled force matrix logically implies a decoupled 
transfer matrix. A comparison of Eq. I56[ Eq. [59] and 
Eq. [5jJ] shows that the force matrix and M s differ only 
in the eigenvalues. Hence all the information that is re- 
quired to compute the decoupling transformation can be 
obtained either from the force matrix or from M s . M c 
can be ignored in the context of decoupling. Even more 
than that: The matrix of eigenvectors E of the matrix M s 
diagonalizes M c , too. We come back to this in Sec. IHI Dl 
and Sec. IIIIF1 after the construction of the matrix of 
eigenvectors E. 

From Eq. [51] one quickly derives that M s is a symplex 
while M c is an antisymplex. As M s and M c share the 
same eigenvectors and since all diagonal matrices, i.e. 73, 
74, 712 and the unit matrix commute, also M s and M c 
commute: 



M s M c - M c M s = . 
From Eq. [15] and Eq. [70] one derives 



M* - m; = 1 , 



which includes that 



= £ C A C + £ S A S , 



(70) 



(71) 



(72) 



in agreement with Eq. 



I. Second Moments and the Envelope Equations 

Besides the position of a beam relative to the design 
orbit, the most important properties of an ion beam are 
collected in the matrix of second moments. We assume 
in the following that the beam is centered, i.e. that the 
first moments are all identically zero. If the state vectors 
ipi( T ) with i = 1 ... n represent a family of n ions, where 
t is the pathlength along the reference orbit, then the 
matrix of second moments a is given by 



1 " 

71 < J 



(73) 



Another possible parameterization is given by a density 
function p(tp, r) = p(qi,Pi,q2,P2, T ) which should be nor- 
malized such that 



J . . . J p{qi,Pi,q2,P2,r)dq 1 dp 1 dq 2 dp 2 



= 1 • (74) 



In this case one writes 



piptp dqi dpi dq 2 dp 2 



(75) 



But independent of the specific practical method of com- 
putation, we assume that the matrix of second moments 
a = ( ip ip T ) is well-defined and has a non- vanishing de- 
terminant. From Eq. [T5]one derives: 



a = (</>V T > + (^ T ) 
= ( F ip ip T ) + (ip i\j t F T ) 

= Fa + a 70 F 70 

We define the S-matrix using xj) = i\) T 70 by 
S = a 70 = ( 4> V> T 7o ) = ( ip i> ) , 



(76) 



(77) 



and obtain from Eq. [76] by multiplication from the right 
with 70: 



FS-SF 



(78) 



In ion beam physics Eq. [75] is called an envelope equa- 
tion, as the second moments define the envelope of an 
ion beam. S is a symplex - as any symmetric matrix 
multiplied by 70: 



T T 

= 7o CT 

= -70 cr 

= 7ocr7o 
= 7oS 70. 



(79) 



Assuming a constant force matrix F, Eq. [78] tells us that 
the second moments are constant, if S and F commute. 
Using the eigenvector-decomposition Eq. [56j of F gives 



= 

= 
= 



FS-SF 

EAE^S-SEAE- 1 
AE _1 SE — E _1 SEA 
AS-SA 



(80) 



Since A is diagonal and S commutes with A, also 

S = E~ 1 SE = D, (81) 

must be diagonal (see Tab. lIVj) . so that the matrix S has 
the form: 



S = EDE 



(82) 



The force matrix and the S-matrix share the same eigen- 
vectors - but will in general have different eigenvalues. 
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J. Matching 

Matching is a concept of circular accelerators, i.e. sys- 
tems with self-feedback, where eigenvalues and -vectors 
are well defined. A beam is matched, if the phase space 
occupied by the ions in the beam fits to the "acceptance" 
of the machine. The practical consequence of mismatch- 
ing is an oscillation or "pumping" of the phase space 
distribution which typically leads to an increase of the 
beam emittance by filamentation 2 ^. 

Wolski formulated the state of matching for the general 
case of periodic motion^. In this case, the restriction is 
no more that F has to be constant, but that F must be 
periodic: F(r + C) — F(r) for a given period C and any 
t. Given the transfer matrix over one turn (or period) is 
M, then the beam is matched, if 

cr = McrM T . (83) 

Using Eq. [T71 one quickly finds in analogy to Eq. [HDJ 

a = -M 0-70 M" 1 70 

0-70 = Mcr7 M _1 

S = M S M 1 (84) 

SM = MS 

= MS-SM 

Following the same arguments as for Eq.|80]onc finds that 
the matrix of second moments is matched, if the transfer 
matrix over one period and the S-matrix share the same 
system of eigenvectors. The general form of the matrix 
S is again given by Eq. [82] where D has a form analogue 
to Eq. [5§] and is given by-ii: 



D = Diag(«£i, -iei,ie 2 , -ie 2 ) , 



(85) 



where the £, are the emittances of the two degrees of 
freedom. 



K. Expectation Values and Scalar Product 

The "expectation values" (ip "f x ip) are related to the 
S-matrix according to: 

$jxip = E ^ k io l i l J^ j 

ijk 

= I E (i> k 10 ii j ¥ + ¥ ii k i ki r) 

ijk 

= \ E (ii j ¥ ^ io l + r ¥ it if) 

ijk 

= h Tr (j2 ill 3 ¥ ^ k 10 + ¥ i> j il k i kl ^j 

= |Tr(7 a ipip + tpip'y x ) 
= |Tr(7 a S + S 7!C ). 

(86) 

The RDM-coefficients of the S-matrix are - apart form 
the sign - the expectation values of the RDMs (see 
Eq.rjl]). 



II. POISSON BRACKETS OF SECOND MOMENTS 
AND THE ELECTROMECHANICAL EQUIVALENCE 

The total time derivative of a function f(p, q, i) is given 
by the Poisson-brackets with the Hamiltonian function 22 : 



df(q,,Pi,t) 
dt 



ill 
dt 



dqi 



df dH 

dqi dpi 



df dH 

dpi dqi 



JLL 

dqo 

df 



70 



V dpi I 



df 
dt 
df 
dt 



( dH \ 
1 d qi » 

dH 
dpi 
dH 
t)q 2 

\ dp 2 J 



(87) 



Vq, P f(p,q,t)tp 

V q , p f(p,q,t)Fip. 



We define the functions fk by the "expectation values" 
according to 

fk(p,q) = \$7kil>- (88) 

fk(p, q) do not explicitely depend on time. Evidently the 
fk vanish for all non-symmetric matrices 70 "fk so that 
there should be exactly n (n + l)/2 = 10 non- vanishing 
functions fk- The gradient V p , g = yields 



(89) 



so that 



^ 7 ViQ 



(90) 



\ (lolkip + ip T lolk) 
\V (707/c + (7o7fe) T ) 
\ V T (70 Ik + Ik 7(f) 

\i? (ioik - ik 70) 

i ip T (70 7^ + 70 707^ 7o) 
§V> (ik+lolllo) ■ 



Since all 7-matrices are either symplices or antisymplices, 
it is obvious, that 



V p,q fk 



ip "f k for k £ [0 ... 9] 
for k G [10 . . . 15] 



(91) 



i.e. only symplices can have a non-vanishing expectation 
value. For all symplices "fk with fee [0 ... 9] the Poisson 
brackets result: 



d_ 

rf7 



iplkip 



ip jk F ip . 



(92) 



On the other hand we find from EQ [18] 



4> T 7o Ik ip + 4> T 7o Ik ip 

ip T F T 70 "fk ip + ip T "folkFtp 

ip T 70 F 70 70 "fk ip + i> T 7o Ik F ip 

-ip F "f k ip + tp "fk F ip 

ip(-f k F -F-f k )ip 



(93) 
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Combining EQ. [92] and [93] yields 
d 

so that 



(ipjkip) = V;2 7fc FV = ^(7feF-F7fc)V, (94) 



= ^2 7fc F^-^(7 fc F-F 7fc )V 

= ^( 7fc F + F 7fc )^ (95) 

According to EQ.JMJ the force matrix F can be written 



(96) 



1=0 



so that EQ. can be written as 



fk — 



d ( ipjk j> \ 
dr \ 2 ) 

9 

J2 Gfe; Fi . 
1=0 



= E V 5 — 

1=0 



! ±^F l 



(97) 



The matrix G^i is composed of the expectation values fi 
and the indices of fi are given by the upper left 10 x 10 
part of the commutator table (Tab. IIV[) . Note that all 
commutators of symplices are either zero or again sym- 
plices. 

Eq. [97] can be reordered such that it has the form of a 
linear transformation of a 10-dimensional vector 



fi - ^2 B ■' ■ 

3=0 



(98) 



The matrix G^z is antisymmetric. The reordering re- 
sults in a quite similar matrix Bjj. Now we introduce 
a normalization that assures a positive sign for positive 
definite second moments. For instance /o is according to 
the definition given by 



jo = ipjoip = i> T ilip 



(99) 



which is negative even though the corresponding sec- 
ond moment is positive. The normalization is done 
by multiplication with 7 ^ = ±1, which equals —1 for 
k e {0,7,8,9}. This can be expressed by the multi- 
plication with a quadratic diagonal matrix g that is an 
extended version of the "metric tensor" g^: 



!J 
9 2 



Diag(-1, 1, 1, 1, 1, 1, 1, -1, -1, -1) 



(100) 



The transformed EQOM are: 



/ = B/ 
gf = (gBg)(gf) 

I = B/, 



(101) 



where the matrix B is explicitely given by 





F 4 


^5 


F 6 


— F 1 


_F 2 


_F 3 












F 9 


-F s 


-F 








F 3 


-F 2 


*5 


-F 9 




"7 




-F 




_F 3 






F 6 


F a 


-F 7 








~ F 


F 2 


-F 1 




-F 1 


F 








F 9 


~ F S 




F 6 




-F 2 




"Q 




-F 9 




F 7 


-F 6 




F 4 


-F 3 






F 


F S 


-F 7 




F 5 


-F 4 








_F 3 


"1 




-F S 


F 5 




F 9 


-F 8 




F 3 




-F ± 


F 6 




-F 4 


-F 9 




F 7 




_F 2 


F l 






F 4 




F s 


-F 7 





Eq. llOll is another way to express the envelope equations 
(Eq. [75] and Eq. ITS]) . The explicite relation between the 
second moments Oij and the expectation values is given 
in App. [El 



A. Symplectic Electrodynamics 

The upper left 4 x 4-block of the matrix B equals the 
electromagnetic field tensor, if we replace the "vector" 
(F 4 , F 5 ,F 6 ) T by the electric field E and (F 7 , F 8 , F 9 ) T by 
the magnetic field B. In this section we investigate the 
equivalence of two-dimensional symplectic flow and rela- 
tivistic electrodynamics. In the next section we show that 
rotations and Lorentz boosts in Minkowski space are for- 
mally identical to a subset of symplectic transformations 
in two-dimensional coupled linear optics. 

If one writes the 4-potential <fr (4-current J, 4- 
momentum P) as a 4-vector using the RDMs 7 o • • ■ 7 3 
according to 

* = ^7o+A B 7i+i4 lf 7a+A.73 = 07o +7^, (103) 
the 4-derivative D as 

D = d t 70 - d x 7 i - d y 72 - d z 73 , (104) 
and the electromagnetic fields E and B as 

F = E xl4 + E yl5 + E z -f 6 + B x -f 7 + B yl8 + B zl9 , (105) 

then the Maxwell equations (MWEQs) can be written 
(remarkably compact) as: 



F = -D* 
DF = 4ttJ, 



(106) 



with the usual choice of units. 

The Lorentz force can also be expressed by RDMs. If 
the 4-momentum is defined by 

P = £ 70 + p x 71 + p y 72 + Pz 73 = £ 7o + PI > (107) 

- where £ is the energy and p the momentum - then the 
Lorentz force equations can be written as 

£ = p =^ Fp - pF )' ( io8 ) 
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where t is the proper time. In the lab frame time dt 
EQ. [TUB] yields (setting c = 1): 



dr 

7 



m " 

I [SE+pxB 
qjvE 

^- [m 7 E + to 7 v x B 
qvE 

q IE + v x B 



(109) 



d£_ 

dr 
dp 

777 

I dt 

I dt 
dE 
dt 
dp 
dt 



which are exactly the Lorentz force equations. 

We use the formal identity of Eg. [751 and Eq. ll08l to in- 
troduce the electromechanical equivalence (EMEQ) 
of 2-dim. symplectic motion and the motion of charged 
particles in electromagnetic fields as described by the 
Lorentz force equations. The statement of the EMEQ 
is, that the transformation properties of the elements of 
the S-matrix with respect to symplectic transformations 
are formally identical to the transformation properties of 
momentum P and electromagnetic field F in Minkowski 
space. The analogy allows to obtain an overview over 
the components of the force matrix F and to establish a 
meaningful nomenclature for the RDM-coefficients: The 
ten elements of the force matrix are associated with the 
energy and momentum of - and with electric and mag- 
netic field components "seen by" - a relativistic charged 
particle: 



£ 70 + Px 71 + Py 72 + Pz 73 

Ex 74 + E y 75 + E z 7 6 
B x 77 + By 7s + B z 79 . 



(110) 



If the force matrix is split into the electrodynamical F e d 
and the mechanical components F TO , the matrices are 
explicitely given by 



f F m 

—E x 

E Z — By 
Ey + B Z 



By 



E.r 

B x 



Ey — B Z 

—B x 

E, r 



+ 



V 


—B x 


-Ey + 


B z E z + 


By 


/ 


~Pz 


£ — Px 





Py 




-£ — px 


Pz 


-Py 










Py 


~Pz 


£ + Px 


V 


Py 





-£ + Px 


Pz 



B x 

-Ey — Bz 
E Z ~ By 

—Ex 

\ 



(111) 

Note that there are ten force terms in symplectic motion 
and only six electromagnetic field components. Corre- 
spondingly there are ten symplectic transformations in 2- 
dim. symplectic flow, but only six transformations known 
in Minkowski space - rotations about and boosts along 
three axis. One also finds that Eq. [TUT] and Eg. 11021 cor- 
respond to the Lorentz force equations only, if F0, Fl, 
F2 and F3 are zero. In this case B is block-diagonal 
- the first block being the upper left 4 x 4-matrix and 
the second block the lower right 6 x 6-matrix. Coupled 



linear motion in two dimensions does not include such 
restrictions and is in this sense a richer theory. 

Nevertheless let us emphasise that the theory of sym- 
plectic motion in two dimensions does not allow scalar 
source terms, i.e. non-vanishing coefficients of 715 in 
the force matrix. If the MWEQs are derived from two- 
dimensional symplectic motion, then the Lorentz gauge 
is included. To make this clearer, we look at the first part 
of EQ. 11061 which includes a scalar expression dt4> + VA: 



'70 
I7o7 



D * = -(d t 70 - 7V) O7o +7-4) 
drf + VA- {d t A + V4 
(V x A) 714707 
(^0 + Vl)l + £7o7' 



(112) 



B 7i47o7 ■ 



where 7 = (7i,72,73) T - Since the force matrix may be 
composed exclusively of symplices, the scalar term has 
to vanish, i.e.: 



d t i> + VA = . 



(113) 



which is the so-called "Lorentz gauge" . In Sec. Ill Fl it 
will be shown, that the so-called duality rotation also 
has properties which are incompatible to the presented 
theory of two-dimensional symplectic motion. 

In the next section we give a survey of symplectic 
transformations using the EMEQ according to Eq. II 101 
This concept allows to obtain a very transparent and 
systematic treatment of two-dimensional coupled linear 
optics. Of course there is some danger of puzzling the 
formally introduced electromagnetic components of the 
EMEQ with the fields of the accelerator due to the simi- 
larity of notation. 



B. A Survey of Symplectic Transformations 

Symplectic transformations have been introduced by 
EQs. [18] and [28j It was shown that the exponentials 
of symplices are symplectic and that there are ten sym- 
plices (in two dimensions) corresponding to ten free pa- 
rameters in symplectic matrices. This indicates a one-to- 
one relation between the symplices and the correspond- 
ing symplectic matrix. In the following we will classify 
the transformations that are based on single 7-matrices. 
Six of these transformations are known as rotations and 
Lorentz boosts in 3-dimensional space, based on the six 
bi- vectors 74 ... 79 , i.e. the "electric and magnetic field 
components" . We call the four remaining transforma- 
tions phase rotation (induced by 70, i.e. by "energy") 
and phase boosts induced by the spatial basis vectors 
71 ... 73 ( "momentum" ) , respectively. All transforma- 
tions are controlled by a continuous parameter a such 
that the inverse transformation is given by the negative 
of this parameter and a vanishing parameter a — yields 
the unit matrix: 



R(0) = 1 
R(a)- 1 = R(-a). 



(114) 
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The parameter a has a specific physical meaning in the 
EMEQ. In case of rotations it is the angle, in case of 
Lorentz boosts it is the inverse hyperbolic tangent of /3 = 
v/c. 

Before discussing the action of the individual 7- 
symplices, we look again at the Hamiltonian in order to 
show exactly why the transformation matrix should be 
symplectic: 

H(p,q) = il> T AiP = -tl> T 7o FiP, (115) 

and insert the equation for a symplectic transformation of 
the forces F' = RF R _1 , then we find (using jp = ip T 70) 

H(p,q) = — f/iR" 1 RF R _1 R-0 

= -(^R^ 1 ) (RFR- 1 ) (Rip) 

= (^R- 1 ) 7 o7o(RFR- 1 )(R^) (116) 

= V' T A'V', 

that the Hamiltonian is invariant under Lorentz transfor- 
mations (LT) and under the assumption that the spinors 
transform according to 

F' = RFR 1 

V/ = R^> 

V/ T = -^ T 7o R- 1 7o = V' T R T (117) 
=> - 7o R- 1 7o=R T 
70 = R T 70 R ■ 

or - in other words: If a transformation matrix R is sym- 
plectic, then the form (and the value) of the Hamiltonian 
and of the EQOM is conserved. Note that if the force ma- 
trix is transformed according to Eq. 11171 then we find for 
the transfer matrix: 



symplex 7^ is given by: 
R 6 = exp(7 b e) 

— (2k){ ^ (2fc+l)! 

k=0 k=0 

00 ok 00 -> i.xi 

= E(-y?) fc ^r + 7fc E(7 b 2 ) fc (f^ 

k=0 k=0 

{1 cos (e) + 76 sin (e) for 7^ = — 1 
1 cosh (e) + 7b sinh (e) for 7^ = 1 
R^ 1 = cxp(-7 b e) 




1 cos (e) — 7;, sin (e) for 7 2 = — 1 
1 cosh (e) — 7b sinh (e) for 7^ = 1 



(120) 

Transformations with 7^ = — 1 arc orthogonal transfor- 
mations, i.e. rotations, while those with 7 2 = 1 are 
boosts. According to Tab. U 7b with b e [0,7,8,9] pro- 
duce rotations and 7b with b e [1 ... 6] boosts. Hence the 
transformed matrices are 

i a = R _1 7 a R 

= R- x 7aR (m) 

= (c- S7b)7a (c + S7b) 

= c 2 7 a - 7b 7 a 7b s 2 + c s (7a 7b - 7b 7a) ) 

where c and s are the (hyperbolic) sine- and cosine- 
functions according to Eq. 11201 The last term on the 
right vanishes if 7b and 7 Q commute. In this case the 
matrix 7 a remains unchanged: 

7a = 7a(c 2 -7b« 2 ) 

J cos 2 (e) + sin 2 (e) = 1 for 7 2 = -1 
1 cosh 2 (e) — sinh 2 (e) = 1 for 7 2 = 1 

(122) 

If the RDMs 7b and 7 a anticommute, one finds: 



M' 



exp(F' s )= £ ^ 



= E 

k=0 



= R E 

Vfe=o 
= RMR 1 



(Fs) k 



(118) 



R 



Symplectic transformation matrices have the form 



R = exp(Fe) 
R 1 = exp(-Fe), 



(119) 



where the matrix F is a symplex. The effect of a basic 



7a = 7a(c 2 +7b- s )+ 2 CS7a7b 

I 7a cos (2 e) + 7a 7b sin (2 e) for j$ = -1 
| 7a cosh (2 e) + 7 a 7b sinh (2 e) for 7 2 = 1 

(123) 

The elementary symplectic transformations (in 2 dimen- 
sions) are summarized in Tab. HU The "frequency dou- 
bling" is an indication that the transformed elements are 
according to their transformation properties second mo- 
ments. This is usually taken into account by using a 
transformation with the "half-angle" . 

Note that even though we describe here the transfor- 
mation of the matrices instead of the RDM-coefficients, 
this is equivalent to the transformation of the RDM- 
coefficients assuming that the RDMs have the same form 
in all coordinate systems. This is only a matter of nota- 
tion. 
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TO 


71 


72 


73 


74 


75 


76 


77 


78 


79 


c = 




coshe 


cosh e 


cosh e 


cosh e 


cosh e 


„ lE 


cos e 




cos e 


s = 


sine 


sinh e 


sinh 6 


sinh e 


sinh e 


sinh e 


sinh e 


Bine 


sint 


sin e 


TO = c 70-- 




+ S74 


+ ST5 


+ S7 6 


- = 71 


— =72 


—=73 








l\ — c Tl ■ ■ 






+ S~ig 


-=78 


-"70 








-=73 


+ = 72 


7 g = c-72 ■ ■ 


-ST5 


-s-yg 




+ S77 




-=70 




+ = 73 




-=71 


T3 = C-Y3.. 


-s-y 6 


+ s-Y8 


-s-Y7 








-S70 


— = 72 


+ = 71 




T4 = C74.. 




+ s-r Q 








+ = 79 


-=7 8 




-=7 6 


75 


= c 75-- 


+ S72 




+ S70 








+ = 77 


+ = 76 




-=74 


7^ = C76-- 








+ "70 


+ = 78 


-S77 




- = 75 


+ = 74 




T7 = c 77-- 








+ S7 2 




-=76 


+ = 75 




-=7g 


+ = 78 


7 R = C T8-- 




+ s 13 




-371 


+ = 76 




-S74 


+ = 79 




-=77 


7q = C T9 ■ ■ 




-s-y2 


+ST1 




- = 75 


+ = 74 




- = 78 


+ = 77 





TABLE II. Table of symplectic transformations in 2 di- 
mensions, a indicates the rows and b the column: "/' a — 
exp(— 7 6 e/2)7 a exp(— 7i,e/2). If j a and j b anticommute, 
then the result is 7„ = c 7^ + s 7 a 75 where c and s are the 
sine- and cosine-function of e, if 7^ = —1, and hyperbolic sine- 
and cosine- function if 7^ = 1, respectively. If j a and 76 com- 



mute, then 7^ 
(Tab. [IVJ) 



7 a . The table follows the commutator-table 



C. Rotations and Lorentz Boosts 

Lorentz boosts are all transformations induced by 7^ 
with b £ [4... 6], i.e. by the electric field terms, while 
rotations are induced by 7& with b 6 [7... 9], i.e. by 
magnetic field terms. Note that both transformations 
are induced by bi- vectors and not by the basic matrices. 
Both - rotations and Lorentz boosts - are well-known 23 
and the matrix coefficients transform exactly as the phys- 
ical quantities associated by the EMEQ. Hence we do not 
need to repeat the formulas here as all relevant properties 
are summarized in Tab. ILT1 

Nevertheless it is instructive to look at the rotation 
matrices with respect to the transformation of the spinor 
tp. The matrix of the spatial rotation about the z-axis, 
induced by 79, is explicitely given by: 



exp (79 77) = 



/ cos ■ 


sin { 

V 0' 





cos ; 


sin { 



— sin ■ 


cos I 






— sin { 


cos I 



/ 



(124) 



Obviously 79 induces a rotation in the plane of q\ and qi- 
The canonical momenta are rotated accordingly. Note 
that the rotation angle in the associated space is e, the 
corresponding angle in the plane of the canonical vari- 
ables is half of that: e/2. 

The rotation about the a;-axis (77) mixes the canonical 
coordinates and momenta q\ and P2 (pi and ^2): 



cos ■ 



exp (77 • 






V - sin § 




cos ■ 
sin -, 





— sin ; 
cos § 



sin I ^ 



cos I 



A phase rotation is a rotation in phase space, for in- 
stance a rotation in the plane of q± and P2 as described 
by Eq. 11251 Another example is given by the rotation 
about the y-axis as induced by 7$ - a combined phase 
rotation in opposite directions: 



exp (78 |) 



cos 1 


sin 


e 
2 










\ 


-sinf 


cos 


c 

2 




















cos 


e 
2 


-sinf 












sin 


e 
2 


cos 1 


/ 



(126) 

In combination with the phase rotation matrix 70 , this 
gives: 

/ 



exp ((70 + 7 8 ) §) 



exp ((70- 7s) §) 



V 

(1 






cose 
— sine 




1 






sine 
cose 






cose 
- sine 



o\ 




1/ 






sine 
cose J 



(127) 



so that both degrees of freedom are here decoupled. Note 
that in this case, the rotation angle e is not divided by 
two. 

We call the remaining symplectic transformations 
phase rotation induced by 70 and phase boosts induced 
by 71, 72 and 73. This is not quite correct, since we saw 
that also 78 induces a phase rotation. Strictly speak- 
ing, there is only one "real" spatial rotation possible in 
two dimensions - and this is induced by 79. All other 
"rotations" are rotations in phase space and are there- 
fore phase rotations. Nevertheless the rotations induced 
by 77 and 78 are spatial rotations in the context of the 
EMEQ. 



D. Phase Rotation 

The symplectic transformation driven by 70 is a "ro- 
tation" between spacelike components and electric field 
components, represented by 707. 70 and the magnetic 
components 77, 78, 79 are unchanged. Expressing this 
using EMEQ one may write: 



E' = E 

p 1 = cos e p — sin e E 

E' = coseE' + sinep 

B' = B 



(128) 



The rotation matrix is given by 
( 

exp (70 § ) = 



(125) 



V 



cos I 
— sin ■ 





sin ■ 
cos ■ 








cos I 
— sin ; 






sin { 
cos ■ 



(129) 
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As listed in Tab. IIII1 the matrices 70, 74, 75, 76 form an 
alternative basis to 70, 7i, 72, 73- The phase rotation pro- 
duced by 70 changes the "mixing" angle of the standard 
basis and this first alternative basis. 



E. Phase Boost 

The transformations induced by 71, 72 and 73 are 
called "phase boosts." These boosts are the exact ana- 
logue of the Lorentz boosts, but the transformation prop- 
erties of E and p are exchanged, as can be seen from 
Tab. |nl While Lorentz boosts are known to have the in- 
variants E B and E 2 — B 2 , a phase boost has the invari- 
ants pB and p 2 — B 2 . Hence the phase boosts transform 
mass and electromagnetic energy into each other as the 
Lorentz boosts are transforming mass and kinetic energy. 
We will come back to this in Sec. IIIII 



To conclude - the duality rotation does not fit into in the 
concept of the EMEQ. 



G. Real Dirac Matrices and Elements in Coupled Linear 
Optics 

In the following we give some examples of force matri- 
ces as used in linear optics of ion beam physic o 20 i 24 . We 
use the transversal coordinates (x,x r ,y,y') as the pairs 
of conjugate variables (q% , p\ , q% , p%) ■ In ion beam physics 
there is the convention to use the path length s along the 
reference orbit as the independent variable. The force 
matrix of a drift (force free motion) is given by 



70 + 76 





( X \ 




(0 


1 





°\ 


X 


d 


x' 
















ds 


y 













1 


y 




\y'J 




u 








0/ 


\y'J 



(133) 



F. The Duality Rotation 

The duality rotation can be expressed by the matrix 
714, but it is not a symplectic transformation and has 
some surprising and puzzling features. Since 714 com- 
mutes with all electric and magnetic field components 
(see Tab. II V[) . the usual form of the transformation can 
not transform the fields, but only the particle properties 
of energy and momentum. Using the nomenclature of 
Eq. 11111 one finds: 



The transversal terms in a solenoid field: 



R 

R 1 

F 

F m 
F'ed 



exp(7i 4 e/2) 

exp(-7i 4 e/2) 

R F e d R 1 = R 1 F e 

RF m R 

RF m R 1 = coseF,, 
RF ed R = coseF ed - 



R 

+ sin e 714 F r 
sine 714 F ed 



(130) 



The last line of Eqs. 11301 is usually referred to as duality 
rotation, explicitely given by: 



E' = cos e E — sin e B 
B' = cos £ B + sin e E 



(131) 



The duality rotation fails to fulfill all reasonable require- 
ments in the context of the EMEQ. First, it is not sym- 
plectic. Second, the transformation is not a similarity 
transformation and finally, it has different forms for the 
mechanical and the electrodynamical components. 

If the MWEQs are deriveable from the EMEQ and 
symplectic flow, then there is no magnetic charge or cur- 
rent density, but the MWEQs are forced to have the 
"classical" form with 



V x E + d t B 



(132) 



d_ 

c/.s 



l x y 
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x'/K 


= K 


y 




V y'/K J 


V 



-1 
-1 
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f x \ 




x'/K 




y 




V y'/K j 



(134) 

where K = 9 with the solenoid field B$, charge q 
and bending radius p. The bending radius depends on 
the momentum p and is given by 



P = 



P 

qB 



This can be written using RDMs as: 

ij) = K (70 - 79)^, 
A radially focussing quadrupole is described by: 



( x ) 




f 





1 









( - \ 


x'/K 


= K 




-1 













x'/K 


y 













1 




y 


\ y'/K J 












1 


0) 




\ y'/K j 



where K 2 = Mr with g = ^ = ^ so that 
again with RDMs - one finds: 



; T r 7o+7i+76 + 78 , 
ip = K ip . 



The axially focusing quad is given by 



ip = K 



7o - 7i + 76 - 78 



(135) 



(136) 



(137) 



written 



(138) 



(139) 



A horizontal bending magnet: 
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X \ 






px 
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-1 















px 


ds 
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py / 
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0) 
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py' I 



(140) 
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expressed by RDMs: 

1 3 7q + 71 + 7 6 + 7s 
P 4 



4>- 



(141) 



Obviously the scale of the momentum coordinates has to 
be adjusted to all individual element types and changes 
even with their excitation (and the beam energy), while 
the scale of the position coordinates is fix. The scaling 
transformations that are required have the following form 



(142) 



which allows to scale coordinates and momenta sepa- 
rately. This matrix - let us call it R - can be written 

as 



f ax > 














x' 


bx' 







6 








ay 










a 





y 


W J 




Ko 








bj 


\y'J 



R 



a + b 



a — b 



'73 : 



(143) 



and it is symplectic, if 
R T 7 oR = f*±*l- 



a + b 
2 

(a + b) 2 



— b \ 

— 73 J 70 

(a-b) ' 



' a + b - 
v 2 



•73 



To 



or 



= a b 70 
ab = l, 

with the equivalent definition 
(* 

R = exp (73 e) = 



(144) 



V 





e e 








e*7 



(145) 



Hence the phase boost in the direction 73 can be under- 
stood as a transformation that changes the relative scale 
of the classical coordinates and momenta. The other 
scaling transformations are induced by the other diag- 
onal RDMs. In the chosen basis these are 74, 712 and 
715 = 1. Besides 73, only 74 is a symplex and is a legiti- 
mate candidate for a symplectic scaling transformation: 



R = exp (74 e) 
= Diag(e~ e ,e e 



(146) 



while 712 and 715 do not generate symplectic transfor- 
mations: 



exp (712 e) = Diag(e £ ,e e ,e £ ,e e ) 
exp (715 e) = Diag(e e ,e e ,e £ ,e e ) 



(147) 



We note that there are symplectic transformations to 
independently change the scales between q% and p\ and 
between qi and P2 , respectively. But there is no symplec- 
tic transformation to change the scales between the di- 
rections {qi,pi) and (92,^2)- Nevertheless, the rescaling 
that is required to compare force matrices of quadrupole, 
solenoids, dipols, etc. to RDMs can be done by symplec- 
tic transformations, namely by R3 = exp (73 e). 



III. DECOUPLING 2-DIM. HARMONIC OSCILLATORS 

Consider a system of two coupled degrees of freedom 
as it is common in many parts of physics. In the general 
case, the force matrix will be time dependent and so will 
be the symplectic transformation matrix R: 



(148) 




RFR 1 J ip = Fip 



Decoupling now means that the transformed force matrix 

F = RR 1 + RFR 1 (149) 

is either a SYSY or block-diagonal. It may still be time- 
dependent. Eq. 11491 can also be written as2£ 

R = F R RF . (150) 
Given that the transformation matrix has the form 

R = exp(-Gr), (151) 
where G is constant, then one obtains 



R = GR 

F = G + RFR 1 



(152) 



If F = const, then it follows from Eq. 11521 and from the 
fact that R and G commute, that 



F = RFR 1 +RFR 1 +RFR 1 = 
= G RFR 1 + RFR 1 + RF G R 1 
F = GF FG. 



(153) 



Obviously, the form of the Lorentz force equation 
fEq. 1108]) is very common in coupled linear optics and can 
be obtained with a few simple assumptions for symplec- 
tic transformation matrices as in Eq. 1 1491 or for symplices 
as in Eq. 11531 or for second moments (i.e. the S-matrix, 
respectively) as in Eq. [761 

In the following we present a straightforward recipe for 
decoupling constant force matrices. It will be shown that 
this method can also be applied to transfer matrices - 
hence can be used to compute the properties of matched 
beams even in cases where the force matrix itsself is not 
constant. 

At the end of this section we scetch some problems 
related to the general case of time-dependent forces. But 
a comprehensive treatment of the time-dependent case is 
beyond the scope of this article. 
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A. Eigenvalues of the force matrix 

The eigenvalues of the force matrix as expressed by 
Eq.[TTTl are given by 

A = Diag(i wi, — i uj\, i 0J2, —i ^2) 
K x = £ 2 + B 2 - E 2 - p 2 
K 2 = -2£ p(E x B) + £ 2 B 2 + E 2 p 2 

- { Ep) 2 - {EBf - (pB) 2 
ui = y/Kl + 2 JK~ 2 

U>2 



= ^K x - 2 ^K~ 2 
Det(F) = K 2 -4K 2 



(154) 



The system is stable, if the eigenvalues are purely imag- 
inary, i.e. if the eigenfrequencies are real. This is the 
case, if the following conditions are fullfilled: 



K 2 > 
Ki > 2 



IK 2 



(155) 



The force matrix can only be a SYSY, if the eigenfre- 
quencies are equal, i.e. if K 2 — 0. This is for instance 
the case, if E = B = 0. If the vectors p, E and B are 
orthogonal to each other such that p = a 2 E x B with 
some constant a, then it follows 



K 2 = -2£p{E x B)+£ 2 B 2 + E 2 p 2 
= -2£pEB + £ 2 B 2 + E 2 p 2 
= (£B-E P ) 2 , 



(156) 



where E = \E\, B = \B\ and p — \p\- In this case the 
force matrix is a SYSY, if £ B = E p. 

We will show in the following that the standard form 
of the block-diagonal force matrix as it can be obtained 
by symplectic transformations, is given by 



/ a 
-P 

\ 

cji = ± \fa~j3 

UJ 2 = ±y/YS- 





7 

0/ 



(157) 



The Hamilton matrix A = —70 corresponding to this 
force matrix is diagonal: 



A = -70 F d 

= Diag(/3, a, S, 7) 

The RDM-composition of is given by 

F d = £ 70 + By 7 8 + p x 71 + E z 7 6 



(158) 



Pr 
E- 



a-\-(3— 7 — 5 
4 

-a+fS+7-6 
4 

a— /3+7 — S 



(159) 



The vectors p, E and B are orthogonal to each other 
so that Eq. 11561 holds, but with a reversed sign for the 
triple product. Obviously one can exchange £ with B y 
and p x with E z without changing the eigenfrequencies. 
This proves that stable systems are possible with B y > £ 
and/or with E z > p x , i.e. with \B\ > £ and /or \E\ > \P\. 

We define a stable system to be regular or massive, 
if all parameters a, /3, 7 and 6 have the same sign (are 
positive), i.e. if the decoupled Hamiltonian function in 
the new coordinates H = tp T A ip is positive definite: 



H = ap 2 1 +l3q z 1 + 1 p z 2 + 6q 2 



-2 



-2 



(160) 



In a regular system the coefficients of the force matrix 
according to Eq. II 571 yield 



£ 2 
£ 2 



> 



> Pi 



(161) 



£ 2 > El 



If one pair, either (a, /?) or (7, ft) is negative, then the 
Hamiltonian function contains a negative kinetic energy 
term and is indefinite. We then call the system irregular 
or magnetic and it is easily shown that in this case 



> £ 2 

> vl 

> El 



(162) 



B. Decoupling a Constant Focussing Channel in two 
Dimensions 

The RDMs and the presented theory of symplectic 
transformations together with the EMEQ facilitates the 
task of decoupling - since we have already the geometri- 
cal and physical understanding that provides the optimal 
strategy to solve the problem. In Sec. II El we discussed 
the meaning of coupling and it turned out that constant 
symplectic force matrices (SYSYs) do not require decou- 
pling. Instead the next derivative yields a decoupled 
second-order equation. In the following we would like 
to discuss decoupling of EQ.[l8]in cases where the force 
matrix is not a SYSY. Then decoupling is the task to 
find a transformation matrix R such that 



RFR 




(163) 



Since products of symplectic transformations are again 
symplectic, we may find the solution in several steps. 

One distinguishes even and odd matrices, where odd 
matrices couple the two degrees of freedom while even 
ones do not*2£. In the chosen basis eight matrices, namely 
7o, 7i, 73, 74 , 76, 78, 712, 7i5 are even and the remaining 
eight matrices are odd. Since the transformation is re- 
quired to be symplectic, the anti-symplices 710 ■ ■ ■ 715 can 
be excluded. There are six even (70, 71, 73, 74, 76 and 
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78) and four odd matrices (72, 75, 77 and 79) remaining. 
This can directly be seen from the force matrix, which is 
given by 



£ 7o + P 7 + E 70 7 + B 714 70 7 



+ By 

• + Py 



+ Py 

z — E 1I 



(164) 



+ Bz 
+ Py 



The odd components are p y , E y , B x and B z correspond- 
ing to the odd matrices 72, 75, 77 and 7g. 

Products of even matrices result even matrices. Hence 
the transformation to block-diagonal form can be done 
by the exclusive use of odd force matrices. Note that the 
odd symplices all anti-commute with each other. Two 
of them (72 and 75) square to +1, the other two to —1. 
Hence they are a basis of the Clifford- algebra CI (2, 2). 

Even matrices may be used for convenience, but exclu- 
sively odd matrices are required to bring a regular force 
matrix into block-diagonal form. The transformations 
have the form of Eqs. 11171 and 11201 where in the latter 
we replace the parameter r with 



R b = exp(7 b e b /2) 
R^ 1 = exp (-7 & e b /2) 
F' = RbFR, 1 



(165) 



After each transformation, the coefficients have to be up- 
dated according to Eqs. ITUl and ITT1 

Then the recipe for regular systems is as follows: 

1. Use the rotation-matrices R7 (about s-axis with 
e 7 = arctan (|r) ) and Rg (about z-axis with e 9 = 
— arctan—) to align the momentum p along the 
y-axis. 

2. Use R5 with £5 = — artanh^- (Lorentz boost) to 
transform into the rest frame. 



B~ 



and 



3. Use again R7 and Rg (with £7 = arctan 

£g = — arctan \Jb^ ) to align the magnetic field 

along the j/-axis. If this is done, we have B x = 
B z =0. 

4. Use Rs (£8 = — arctan ^J 2 ^) to rotate about y- 

axis and make E x vanish. Note that this step is 
not required to obtain block-diagonal form. But it 
is required for a complete diagonalization. 

5. Use R2 (£2 = artanh^-^-^ ) to make the parallel 
electric field component E y vanish. 

All but the last step should be clear. In the rest frame 
we have p = 0, so that rotation matrices can be applied 
without restrictions. The (odd) 72-term (i.e. p y ) will 
remain zero. The alignment of B along the y-axis brings 
the odd terms B x and B z to zero, the rotation about the 



y-axis yields E x 

75, 



0. The only remaining odd term is 



i.e. E y : 



F (3) = m/yo + E 75 + e z j 6 + By 7s 



By + E z 



-By + E z + , 




(166) 

where we wrote m instead of £ as the 70 coefficient to 
make clear that the transformation into the rest frame 
has already been done. The last step is required to 
achieve the block-diagonal form of Eq. 11571 where 

F< 4 ) = rh-f + E zle + B y -f 8 



a 

7 
S 



By- 
By- 

-B, 



E z + m 
E z + m 
+ E,+m 



(167) 



— — B,. — E, 



m 



where rh is given by 



to 2 — £" 2 . 



(168) 



This shows that phase boosts are transformations which 
convert electromagnetic energy into mass and vice versa 
in the same sense as Lorentz boosts convert mass into 
kinetic energy or vice versa. Note that we could also ex- 
change the order of the last two steps, i.e. of Rs and R2. 
The result would be the same since 78 and 72 commute. 

If the force matrix is irregular (but stable), then the 
energy term might be too small to perform the Lorentz 
boost into the rest frame as suggested by the recipe. In 
this case, the magnetic components provide the focusing 
strength that stabilize the system (see Eg. I154j) . 

In this case we have to prepare the force matrix in the 
following way: 

1. In the irregular case, we have to consider the pos- 
sibility that p 2 > £ 2 . In this case the first step is a 



phase rotation Rq with eq = j arctan ^-J^-S^ 



to 



minimize p 1 . 



2. Align B along the y-axis. 

3. Rotate about y-axis to maximize p z (i.e. to make 

Vx = 0). 

4. Apply phase boost using 71 by £1 = \ artanh(^). 
If the energy is still to low, perform the following steps: 

5. Align (again) B along the y-axis. 

6. Rotate about y-axis to maximize E z (i.e. to make 
E x - 0). 

7. Apply Lorentz boost using 74 by £4 = 
iartanh(fj). 
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Now the energy £ should fulfill £ 2 > p 2 , so that the usual 
recipe for regular matrices can be applied. 

An inspection of Tab. IIIII may help to explain these 
transformations. Magnetic force matrices are irregular 
with respect to the definition of the state vector, i.e. the 
choice of 70 as explained above with Eq. [1] and Eq. [5J If 
we inspect Tab. IIIII then we can see that in basis system 
9), where 7§ represents the energy and 73 represents p x , 
7! is the driver of a Lorentz boost. Interpreted in the 
usual system this means that it is possible to minimize 
p z using Ri, if B 2 > p 2 . The second suggested trans- 
formation corresponds to the same basis and a Lorentz 
boost along y (accordingly in basis system No. 9, this 
is driven by 74). The fact that we have to look at a 
different basis system to find the appropriate symplectic 
transformation, legitimizes us to claim that such systems 
are irregular. The examples given in Sec. Ill Gl have 70- 
cocfficicnts which are at least as strong as any other com- 
ponent. This fact seems to support the assumption that 
only regular systems exist in linear coupled optics. This 
assumption is wrong and in the following section we give 
an example for an irregular system. 



C. Example for an Irregular System 

In an accompanying paper we describe a simplified and 
idealized cyclotron model with space charge, which is an 
example for an irregular system^. The force matrix can 
easily be transformed into the form of Eq. 11671 Using 
ip = (x,x',l,5) T as the canonical coordinates, where x 
and I are the horizontal and longitudinal position of the 



ion, x 



dx 
ds 



is the horizontal angle and 5 



po 



is the 



relative deviation of the momentum p from the average 
momentum po- The (constant) force matrix F is given 
by 



V 



100 

-k x + K x h 
-h -4 



\ 







(169) 



K zl 2 / 



The RDM-coefficients are 

£ = 

P x = 

Py = 

E x — 

Ey = 

E z = 

B„ = 




1 2 K Z 
1 2 K Z 



(171) 



1 Z K 2 



The system is resonably simple as there are only two ele- 
ments that are not block-diagonal. This allows to guess a 
transformation matrix instead of using the above recipe: 



R = exp((r + l/r)|7 2 + (r-l/r)f 77 ) 
= cosh (s) 1 + ((r + l/r) 72 + (r - l/r) 

Where the follwoing abbreviations are used: 



(172) 



A — h 



B 

cosh (s) 
sinh (s) = — 



u> 2 +K z 



B-A 
~A 



(173) 



B-A 



The transformed force matrix is then 



RFR 











1 2 K Z 

_ AK Z +B (fcg-iCc)- 





2 h 



(174) 



B-A 

B K z +A(ky-K x )-2h 
(B-A)^ K z 



Since 7 2 K z > and T > is given, the transformed force 
matrix is irregular and the transformed Hamiltonian H 
has the form: 



H 



x' 2 



+ /3x 2 - TS* -i*K x F 



(175) 



where 7 is the relativistic factor, h = l/r is the inverse 
bending radius of the magnetic field and k x = h 2 (1+n) = 
h 2 7 2 is the horizontal restoring force, n — -g ^ is the 
field index. K x and K z are the horizontal and axial space 
charge forced, respectively. The eigenfrequencies of the 
force matrix are 



It is also an example for a classical system, that can not 
be diagonalized by the method of "symplectic rotation" 
of Teng and Edwards^ 2 -. 



D. Diagonalization 



= kx-Ky-Kz Though the force matrix has been decoupled, there is 



ui 



b = K z {K x + h 2 1 2 



. motivation to continue the process of diagonalization by 

'' (170) a f ew m ore transformations, since then the product of 

^ = V a + v 1 °- 2 ~ b the transformation matrices equals the matrix of eigen- 

\J a — \J a 2 — b . vectors and the resulting diagonal matrix contains the 
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eigenvalues. According to Eg . l56l l57l l82l and Eg. 1551 tne 
matched beam-ellipse can then be directly computed for 
any combination of the emittances E\ and 

The intention is to bring F into the form of EQ. [56] 
Since we expect imaginary eigenvalues, the transforma- 
tion matrices can not be real-symplectic. The first "direc- 
tion" that we use is i 75 , which can formally be identified 
with a Lorentz boost with an infinite imaginary relative 
"velocity" , so that the "angle" is given by 



so that 

F0O = P 8 F( 6 )P^ = A 



a = 
b = 




= Jl- 



1 - 





i (By + m) c 





El 



(By 



El 
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i (By -fh)b 

i (By - m) b 



4a/3 
4 (57 

(3+tF 



lim artanh(i /3)/2 

0— >oo 

lim I arctan (/J) 
i - 

1 . 
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(181) 

In terms of RDMs the diagonalized force matrix is: 

F (7) = -i[a(B y + m)-b(By-m)}~/ 3 

- Ua{B v + m) + b(B y -m)]~ii (182) 



The matrix is explicitely given by 

P 6 = exp(i7 6 f) 

/ 1 -i \ 
-i 1 
ooi-?: 

\ -% 1 J 



1 

V2 



(177) 



This matrix has some remarkable properties: It is uni- 
tary, symmetric and symplectic. The transformed force 
matrix is 



p(6) 



PeF^Pg 1 






i (By - , 



i(By 



(178) 

And the final transformation includes two RDMs, 70 and 
78 . This combination allows for a relative phase rotation 
between both degrees of freedom as in EQ. I127( but now 
with imaginary angles: 



exp [1 70 
/ Cl 
-isi 




ei-e 2 
2 

i si 
ci 

C 2 



\ is 2 



-^78^) 
\ 


-i S2 
C2 



(179) 



Note that the last (double-) transformations is not nec- 
essary, if E z is already zero. The eigenvectors are now 
given by the product of all transformation matrices (in 
their order). Since Rs is also symplectic, it turned out 
that the matrix of eigenvectors E as a product of (com- 
plex) symplectic matrices is (complex) symplectic. 



E. Summary: Decoupling and Diagonalization 

The force matrix F of stable systems was shown to have 
the form of Eq. [SB] where the diagonal matrix A has the 
form of Eq.[5j|l It turned out that there is an intermediate 
stage in the process which corresponds to a diagonalized 
Hamiltonian function. The similarity transformation R 
to reach this intermediate stage is symplectic and with 
Eq. H7] and Eq. HH it follows that 



A d 
F d 

7o A d 



Diag(/3,a,(5,7) 
RFR 1 
R70 AR 1 
-R70 A7 R T 7o 
70R70 A7 R T 7o . 



(183) 



This equation can be reversed by multiplication with R 
(R T ) from the right (left), respectively: 



A = R A rf R . 



(184) 



where 



Note that R is in general not orthogonal, ie. R T R^l. 
Hence a symplectic diagonalization can usually not be 
replaced by an orthogonal transformation. Nevertheless 
the symplectic decoupling diagonalizes the Hamiltonian. 



Ci 
Si 

£2 



cosh Si 
sinh Si 

i A rtanh(-^- 
1 Artanh( 7r ^- 



(180) 



F. The RDM-Coefficients of the Transfer Matrix and the 
Tunes 

In the following we use the knowledge of the matrix of 
eigenvectors E to further investigate the general structure 
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of the transfer matrix using Eq. [66] and Eq. [69] The 
"amplitudes" of M s arc limited by the sine-function, but 
the structure is identical to the structure of the force 
matrix, i.e. M s is composed of symplices. The structural 
difference between transfer and force matrix is given by 
M c . 

Note that the decoupling algorithm of Teng and Ed- 
wards for twodimensional systems is applied to the trans- 
fer matrix instead of the force matrix. But from Eq. 1661 
and Eq. IH3] its it clear that the computation of the eigen- 
vectors does not require to include M c . The decoupling 
can therefore be done with the force matrix F, - which is 
favorable in case of constant forces - or with M s . The lat- 
ter enables the treatment of periodic systems with vary- 
ing force matrix F. If the decoupling algorithm is applied 
to M s , then one obtains by comparison of Eq. Eq. rjj 
and Eq. [1671 



A. 



\ Tr(M) 

a(a+P)+b(-y+S) 

q(q+/3)-b(7+<5) 
4 



(185) 



The last missing value to determine the coupled tunes is 
A c . We will argue in the following, that the decoupling 
of M s will casually diagonalize M c and hence yields A c . 

It is noteworthy that 712 is the only even matrix that 
anticommutes with all odd symplices (see Tab. IVT). This 
means that 712 plus the odd symplices 72, 75 and 77 
build an alternative basis of C7(3, 1) - in which 79 is the 
pseudoscalar (see Tab. IIII[) . 

In the following we use the diagonalization steps and 
reconstruct the components of M c . The scalar part can 
be ignored as it can not be transformed. The interesting 
part therefore is E712 E _1 . The matrix E is the reverse 
product of all transformation matrices and hence we may 
analyze the structure of M c by applying all steps of the 
diagonalization in reversed order. 

As can be seen in Tab. IIV1 the matrix 712 commutes 
with 70, 76 and 7s, so that the reversed transformations 
of Eq. 11781 and Eq. 11811 have no effect on 712. Hence 
we can skip the diagonalization steps and go back to the 
transformations of Sec. lIII Bl The general structure of the 
resulting terms can be seen from Tab. IIVI All transfor- 
mations are symplectic, i.e. are launched by symplices. 
All non-vanishing commutators of RDMs in the range 
7io • • ■ 714 with symplices yield matrices in the range 
7io ■ • • 7i4- Therefore the symplectic (back-) transforma- 
tion of 712 exclusively yield axial- vector and pseudoscalar 
components. In the following we assume that the one- 
turn-transfer matrix is known and "off resonance" , i.e. 
non of the coefficients of Eq. [S3] are zero. 

Then the first step of the back-transformation is 



X« = R^7i2R2 

= cosh e 2 712 + sinhe 2 7i4 



(186) 



The second step - the rotation about the y-axis using 73 
again has no effect, since both, 712 and 714 commute with 



78- Hence the next steps are the rotations about z- and 
x-axis. Both rotations commute with 714 - so that the 
second term in Eq. 1 1861 remains unchanged. But the first 
term is an (axial) vector term which rotates analogue to 
the vector term 72. Hence the rotations create coefficients 
of 711 and 713. The next step is the inverse lorentz boost 
in y-direction with 75. Since 75 commutes with 711 and 
713, only the 712-term transforms and launches a 7x0- 
component. 

The remaining last two transformations are again ro- 
tations about the x- and z-axis. Rotations will neither 
change the pseudoscalar coefficient nor the time-like ax- 
ial vector component 710. Hence the pseudoscalar coef- 
ficient of the transfer matrix can be used to determine, 
if a phase boost is required for decoupling. According to 
the decoupling-recipe this is the case, if (and only if) the 
EMEQ-components of the force matrix fulfill EB = 0. 

Hence one finds: 



— sinh £2 A c = 7x4 ' M 



X = M c — ~ Tr(M c ) 
= — A c E 7x2 E _1 



(187) 



If the decoupling transformation has been computed, 
then 

A c = - 7l2 E- 1 (m c - ~ Tr(M c ) J E . (188) 
Hence we may now solve Eq. [53] for the tunes: 



tan (lot) 
tan (Awr) 



E, 
£<■• 
A, 



Ag 



S 8 



(189) 



If t is the repetition time (or length of the design orbit) 
of the accelerator, then the tune Qi equals Qi = 

Note that the anticommutators 7> ' 7,y "^ 7 " 7> ' for /j,, v 6 
[10 . . . 14] form a diagonal 5 x 5-matrix and hence are the 
basis of a five-dimensional phase space (see Tab. IV]) . 



G. Transformation to Equal Frequencies 

Consider now a time dependent symplectic transfor- 
mation as derived in Eq. 11521 with 



G = i lo 74 , 
with the corresponding matrices 
U = e l " 74T 

U — 1 = g~ '<*>74T 

= Diag(e' i " T ,e~' iWT ,e iwr ,e- il ' jr ) 



(190) 



(191) 



which does not affect the diagonalizcd force matrix F: 
UFU _1 =F, (192) 
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but results in a frequency shift such that 



F' 



9-2 



F + i uj 74 

\ 
-ifli 
iVL 2 
\ — z 2 y 

(m + By ) a — uj 
(rh — By) b + to 



(193) 



One may now chose lo such that Qi = fig = fi, that is: 

(rh + By) a — uj = (rh — B y )b + uj 

2uj = (rh + By) a - (rh - B v )b (194) 

rh (a—b\+Bv (a-\-b) 
UJ = 1 2 



then we obtain the result that 

F = -i ft 73 . 
It is also possible to chose 

G = i uj 73 , 
resulting in a frequency shift such that 



(195) 
(196) 



F' = F 



■ i uj 73 



-i ^74 

ft (a+b)4 
fi = (m + By)a — uj . 



m (a+b)+i3 y (a—b) 

U - 2 



(197) 



The general time dependent case is more involved. If 
we intend to use Rf, in order to let the corresponding 
odd component of the force matrix vanish, then Eq. 11991 
implies that we have to take the new (odd) term 75 i into 
account. 

Following the recipe from above, the first step is the 
alignment of the momentum along the y-axis using R7 
and R 9 , respectively. This introduces new (time depen- 
dent) terms in the coefficients of 77 and 79 - which is no 
problem at this stage. The second step is a boost along 
the y-axis using 75 with the intention to get rid of p y . 
Even this is not a problem, since 75 corresponds to E y , 
which has not yet been considered. Hence it is always 
possible to transform into the "rest frame" at the cost of 
additional "field" terms in B x , B z and E y . 

Nevertheless the next steps include difficulties, since 
the alignment of B along the y-axis induces again addi- 
tional terms in B x , B z . If the generators of the transfor- 
mation and the terms to act on are identical, the straight- 
forward approach fails. The rotation about the x-axis: 



arctan 

B z B v — B v B z 



(A) 



B, 



B x 



More general one has 



H. Time Dependent Forces 



Bl+Bl 

I B z By— By B z 



arctan (y\ 

xy-y x 

Y 2 +X 2 

artanh(Y") 

XY-YX 
Y 2 -X 2 ■ 
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(204) 



The treatment of time dependent symplectic transfor- 
mations goes beyond the scope of this article. Neverthe- 
less we would like to make a few comments. 

If the parameter e of Eq. 11201 is time dependent, it 
follows that 

R 6 R^ 1 =7 b e. (198) 

Eq. 11491 can then be written as: 

F = 7fc e + R h FR^ 1 . (199) 

Consider for instance the case of a force matrix F with 
the simple periodic form 

F = F (l+acos(wr)). (200) 

Then the use of a symplectic transformation matrix 

R = cxp (—a/w F sin(w r)) , (201) 

results in a constant force matrix: 

tp = (— aF cos(wt) + F + F a cos (lot)) tJj / 9n9 \ 



IV. SUMMARY 

We introduced the real Dirac matrices and the concept 
of the symplex in two-dimensional coupled linear optics, 
i.e. classical Hamiltonian mechanics. Since these ma- 
trices form a basis for all real-valued 4 x 4-matrices, we 
can claim that our survey of symplectic transformations 
for such systems is complete and final. We have shown 
that a subset of these transformations is isomorphic to 
rotations and Lorentz boost as applied to Dirac spinors 
in quantum electrodynamics. We used this isomorphism 
to describe a general and straightforward decoupling al- 
gorithm for the case of constant forces and for transfer 
matrices in the context of periodic motion. The algo- 
rithm was tested numerically using random force matri- 
ces of stable systems and it never failed in case of stable 
systems. 

There are six possible choices for the symplectic unit 
matrix 70 as listed in Tab. IIIII corresponding to six an- 
tisymmetric RDMs. These choices are nothing but per- 
mutations of the dynamical variables in the state vector, 
and hence they do not differ in their physical content. For 
each possible 70 there are two choices for 71, 72 and 73, 
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which can be transformed into each other by the phase 
rotation. Hence the choice of the system of RDMs is 
unique as soon as the order of the canonical coordinates 
and momenta is fixed. 

We clarified the structure of the transfer matrix in 
two-dimensional coupled linear optics which we believe is 
the first universally valid generalization of the Courant- 
Snyder theory. 

V. DISCUSSION AND OUTLOOK 

Since the Dirac matrices are usually known in the con- 
text of relativistic quantum electrodynamics (QED), it 
is sort of inevitable to contemplate the relationship be- 
tween QED and classical mechanics. It is known that the 
eigenvalues of stable systems lie on the unit circle in the 
complex planed, so that the treatment of eigenvectors - 
which are essential to the theory - requires the introduc- 
tion of imaginary numbers. David Hestenes gave a geo- 
metrical interpretation of the unit imaginary in quantum 
mechanic o 27 i 28 . In context of two-dimensional symplectic 
flow, the unit imaginary can be interpreted statistically, 
as the complex eigenvectors are able to replace an statis- 
tical ensemble of real- valued trajectories. This is to say, 
that the matrix a — ipip with a real-valued "spinor" ip 
has a vanishing determinant, i.e. zero emittance - un- 
less we introduce some kind of averaging - either over 
a statistical ensemble or over time. The introduction of 
complex- valued eigenvectors circumvents this problem as 
well-known in the Cour ant- Snyder theory. 

There is a theorem that all symplectic systems of 
the same dimension are isomorphic 18 . If the reverse 
is also true then the existence of the EMEQ suggests 
that Minkowski space-time is "embedded" in a four- 
dimensional phase space and that special relativity can 
be "deduced" from classical symplectic motion. This de- 
duction implies that the components of the relativistic 
momentum are not proper canonical variables in the clas- 
sical sense, but instead energy and momentum appear 
to be second moments of the dynamical variables of a 
two-dimensional symplectic system. In fact, the invari- 
ance of mass is then formally and physically identical to 
the invariance of the phase space volume (emittance) in 
symplectic flow. This formal analogy seems to find its 
correspondence in the relation 



E 
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There are still many open questions as for instance: 
why are there six held components in electromagnetic 
theory instead of ten as we should expect from two- 
dimensional symplectic motion? The answer could be 
related to the dimensionality of the system under con- 
sideration. The phase-space dimension n of a classical 
system and the number v of "valid" symplices are re- 
lated by Eq . |2U1 and can be resolved for n: 



2v 



1 

4 



If we insert v = 6, then we find n = 3. Classically the 
phase space dimension should be an even number, since 
Hamilton mechanics is based on pairs (qi,Pi)- But a gen- 
eralization of Hamilton mechanics for odd phase space di- 
mensions has already been presented by Nambu utilizing 
several distinct Hamilton functions 2 ^. We believe that 
the description of odd-dimensional phase spaces can also 
be done if canonical pairs (qi,Pi) are used - if the dimen- 
sionality of a systems is defined by the number of avail- 
able force components, i.e. the number of symplices. Rel- 
ativistic electromagnetic theory knows six symplices (and 
hence six symplectic transformations), two-dimensional 
coupled optics has ten, in App. [A] we argue that com- 
plexification of the presented theory leads to a 15 valid 
symplices corresponding to a 5-dimensional phase space. 
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Appendix A: Complexification 

If we consider complex spinors then it is quickly shown 
that the Hamilton equations are: 

L _ q l +i P i 



n - J— { ML _ ML 

» y/2 dip* 

j) - l (ml _i_ ml) 

Pi — ^2 V W>< J 
OH 



(Al) 



4 # = -H 



Given the Hamiltonian function is expressed by 

H = tfAil> = Y,1>i-Aij1>j, (A2) 



with the matrix A then one obtains: 

= -J2i>iA ik 

i 

With implicit summation we write 



(A3) 



i ip = Aip 
itpi = -ipi A 



(A4) 
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If we take the adjunct of the second equation 

iiP = {A T )*ij = A i ^, (A5) 
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it becomes obvious that the matrix A must be self-adjoint 
or hermitian. If we switch back to the (qi,Pi) approach 
with A = A + i B, A = A T , B T = -B 



H 
H 
H 



ip^ Aip 

\{q-ip) T (A + iB) (q + ip) 
I (q T Aq+p T Ap-2q T Bp) 



(A6) 



then we find additional restrictions compared to a "classi- 
cal" system: The matrix A appears for both - coordinates 
and momenta. The complex notation of a 4-dimensional 
spinor is therefore not suitable for the general description 
of 4-dimensional classical harmonic oscillators. 

The symmetric matrix A in the real- valued 2-dim. de- 
scription has been replaced by an hermitian matrix in 
the complex notation. It can be shown that the gener- 
alization of a symplectic transformation matrix U then 
is 



To 



= U 7o U t 



(A7) 



There are 15 complex Dirac matrices with zero trace that 
obey Eq. IA7I so that the phase space is - according to 
Eq. 12061 5-dimensional. 



Appendix B: The 7-Matrices 

To complete the list of the real 7-matrices used 
throughout this paper: 
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TABLE III. List of the indices of equivalent canonical basis 
sets of 7-matrices and the corresponding re-interpretations 
of the matrices relative to the choice of the basis. Negative 
indices indicate the negative matrix with the index taken as 
absolute value. (Vo, V) are the elements corresponding to the 
so-called "axial vector" , "PS" is the pseudo-scalar. 



Appendix C: Arbitrary Hamiltonians 

The concept presented above is far more general, as 
any arbitrary Hamiltonian H(*£) may be written as a 
Taylor serie up to second order: 

= h + J2 ^ ** + \ Aik **+•••• ( C1 ) 

k i,k 

From Eq. ED and Eq. M one finds the following EQOM: 
* = 70 £ + 70 A * . (C2) 

Leach formulated a time canonical transformation to 
transform a time-dependent Hamiltonian of the form 
given by EQ. ICll to a different (simpler) form, which 
we will summarize in the following. Given a symplectic 
time-dependent transformation matrix R, one writes for 
the spinor ip in the transformed system: 



ip = R * + ip . 
The new Hamiltonian H should have the form 
H = H +Hj + ^ t Aip. 

and the new EQOM are 

ip = 70 e + 7o A ip , 
or - using EQ. El and EQ. [C2j 



(C3) 



(C4) 



(C5) 



ip = R* + R* + ip = R * + R (70 e + 70 A *) + ip . 

(C6) 

Combining EQ. Eland EQ. ES yields: 



ipo = 70 Ai/j - R7 £ + 7o£ 
R = 70 A R — R 70 A , 



(C7) 



where the "coordinate-free and coordinate-dependent 
parts have been separated"—. 
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We are especially interested in a transformation which 
keeps the forces 70 A = 70 A but sets e — in the new 
system. One obtains with F = 70 A: 



ip = Ftp 

ipo = Ftp - R,7 i 
R = FR RF, 



(C8) 



If one considers the (rather trivial) situation that e = 0, 
then the conclusion is that any spinor ip can formally be 
replaced by ipo + R ip - or more generally: Given a set of 
vectors ipk where each vector ipk fulfills the EQOM 

Tpk=~Fipk (C9) 
then one finds that a superposition of the form 

ip = J2 Rfe i> k (cio) 

k 

where the symplectic matrices R& obey the EQOM 

R*; = F Rfc — R*: F , (Gil) 

also fulfills the EQOM: 

ip = R *; 4>k + R fe ^fc 

= E^F^ + E (FR fe -R fe F) ^ fc 

fc fc (C12) 

= F0 

Appendix D: Multiplication Tables of the RDMs 

Tab. IIVI lists the commutators and Tab. |V] the anti- 
commutators of all RDMs. The sum of both is the mul- 
tiplication table of the group. 

Appendix E: Second Moments and Expectation Values 

In order to obtain the direct relation between the sec- 
ond moments cry according to Eq. l73l and the expectation 
values fk, we established two new vectors. The first is a 
vector of second moments: 



3 = (fll, 022, 033, 044, 012, 034, 013, 024, 023, 014) 



(El) 



and the second is the vector / of the form 

x = (/o, fi,fe, /s, /3, A, /2, JV, /5, fn) T ■ (E2) 
Ordered in this form, the equation 

a = Tx (E3) 



TABLE IV. Commutator table of the group of 7-matrices: 
^ _ 7a 7b -76 7a ^rjjgj-g a j g ^jjg i nc j ex Q f the row, b the index of 

the column and c the table entry. The only matrix commuting 
with all other matrices, is the unit matrix 715 = 1 



12 3 4 



9 10 11 12 13 14 1 



TABLE V. Anticommutator-table of the group of 7-matrices: 
_ 7a 7b+7 7a wnere a j s the index of the row, b the index 
of the column and c the table entry. 



holds for a block-diagonal matrix T: 

/ 1 -1 1 1 

II- 11 

III- 1 
1 -1 -1 -1 



1 1 

1 -1 



1 1 
-1 1 



(E4) 



-1 -1 
-1 1 / 



Note that the ordering is chosen such that the first six 
(last four) elements of x correspond to expectation values 
of even (odd) RDMs. The inverse T _1 has the same 
block-structure. 
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